Paramagnetic effect in YBa 2 Cu 3 07_ :r grain boundary junctions. 
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A detailed investigation of the magnetic response of YI^CusOz grain boundary Josephson junc- 
tions has been carried out using both radio-frequency measurements and Scanning SQUID Mi- 
croscopy. In a nominally zero-field-cooled regime we observed a paramagnetic response at low 
external fields for 45° asymmetric grain boundaries. We argue that the observed phenomenology 
results from the d-wave order parameter symmetry and depends on Andreev bound states. 

PACS numbers: 74.50.+r, 74.40.+k, 85.25.Cp, 73.23.Hk 



Superconductors below the transition temperature T c 
usually expel an external magnetic field. This phe- 
nomenon, known as the Meissner effect, leads to diamag- 
netism. It was therefore quite unexpected that a para- 
magnetic Meissner effect (PME) was observed in a field 
cooled regime for ceramic P^S^CaC^Os 1,2 (BiSCCO). 
For an explanation, it was proposed that below T c there 
are randomly distributed spontaneous orbital currents in 
the samplei^ These currents can be oriented by an ex- 
ternal magnetic field, providing a paramagnetic signal. 
Sigrist and Rice 4 pointed out that the spontaneous or- 
bital currents can be caused by a d x 2_ y 2-w&ve symmetry 
of the superconducting state. Indeed, the d-w&ve scenario 
predicts the existence of Josephson junctions in which the 
Josephson energy reaches a minimum for a phase differ- 
ence across the junction at <j> = 7r (7r-contact). In ceramic 
high-T c superconductors (HTS) the grains can form loops 
which contain odd numbers of 7r-contacts, so called frus- 
trated loops^ In such a system the energy is minimized 
by a configuration with a spontaneous current flowing 
in the loop£ In particular, for a 7r-contact with a "con- 
ventional" current-phase relationship I = I c sin(0 + it) 
inserted in a loop, the gain in Josephson energy exceeds 
the loss of magnetic energy when the parameter j3 satis- 
fies f3 = 2irLI c /$o > 1, where $o is the flux quantum, I c 
is the critical current and L is a loop inductance. How- 
ever, a paramagnetic signal in the field-cooled regime has 
also been observed for conventional superconductors!^ 

In order to distinguish between different origins of 
paramagnetism Rice and Sigrist proposed to detect 
spontaneous orbital currents by a SQUID microscope for 
zero-field cooled samples. In a granular BiSCCO sam- 
ple exhibiting a paramagnetic signal, spontaneous mag- 
netization has indeed been observed^ The PME and the 
presence of spontaneous currents therefore represent two 
of the main features induced by the d-wave order param- 
eter symmetry in HTS. 

A study of these phenomena on a more controlled sys- 
tem, such as a grain boundary (GB) line, could shed 
additional light on the relation between these effects and 



their influence on the properties of the grain boundary 
Josephson junctions (GBJJ). 

In this paper, for the first time, we give evidence 
of a paramagnetic behavior for zero-field cooled asym- 
metric 45° YBazCugOz (YBCO) GBJJ. The paramag- 
netic signal, that has been previously observed in two- 
dimensional systems, has been in this case detected in a 
single GB line and related to basic mechanisms in HTS 
JJs. Apart from conventional transport properties and 
Scanning SQUID Microscopy (SSM) measurements, we 
exploited radio-frequency measurements, which are very 
sensitive and provide a direct test for paramagnetic sig- 
nals as demonstrated below. This is a further manifesta- 
tion of the e?-wave nature of the order parameter, but it 
is also relevant for the understanding of transport across 
GBJJs and in particular on the incidence of Andreev 
bound states. 

The idea of the rf measurements is the following. The 
sample of interest is inductively coupled to a high-quality 
(quality factor Q w 300) parallel resonant circuit of in- 
ductance Lt, capacitance Ct- Experimentally this is 
realized by flip chip configuration - the sample is placed 
on the top of a small solenoid coil perpendicular to its 
axis. For such an arrangement, the effective impedance 
of the tank circuit coupled to the sample is a function 
of the external magnetic field H e applied to the sample. 
This field consists of a dc and an rf component as in- 
duced by a current Id c (in fact of very low frequency) 
and an rf oscillating current I r f in the tank coil, hence 
H e = Hdc + H r f . In order to avoid the nucleation of rf 
"vortices" and other nonlinear effects, the amplitude of 
H r f is small, so that H r f <C H c \,Hd c where H c \ is the 
first critical field, therefore H e = Hd c - If L e ff(Hd c ) and 
R e ff(Hdc) are the effective inductance and the effective 
resistance of the tank circuit-sample system respectively, 
the phase angle a between the drive current I r f and the 
tank voltage U is given by 

tana = r — 7TF~) ' ( ~7^~ ~ uL ef f{H dc )) ■ (1) 

Hef f (Hdc) \UJCt ) 
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It follows from Eq. ^ f° r Hdc 
frequency luq 



and at the resonance 



1 / y/L e ff(0)CT that the parameter a is 
zero. Therefore, by monitoring a as a function of Hd c 
at the frequency u>o, the L e ff(Hd c ) dependence can be 
obtained (note that R e ff(Hd c ) dependence is controlled 
independently by measuring the quality factor Q of the 
tank circuit-sample system). If R e ff — Rt does not 
depend on Hd c then tan a can be written as: 
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where k is the coupling coefficient between the tank coil 
and the sample, and \m is ac magnetic susceptibility of 
the sample. In the experiments we are obviously inter- 
ested in the regime Xm ^ const, the only configuration 
able to provide significant information. Let us analyze 
Eq. [21 near H e = 0. When the supercurrent induced by 
an externally applied field is diamagnetic, a change of 
Xm < will result in a local maximum in a(H e ). Sim- 
ilarly, a local minimum of the u(H e ) curve indicates a 
paramagnetic response (Xm > 0). 

The simplest system which exhibits a similar behavior 
is the well-known rf SQUID. The sensor of this device is 
a Joscphson junction inserted in a superconducting loop. 
If the inductance L of the loop is relatively small (so that 
[3 < 1), the a(Hdc) dependence has a local maximum at 
Hdc = (see for example Ref. Ull). Due to the induced 
current in the sensor, the magnetic flux inside the loop 
satisfies < <i> e , where $ e is an applied external flux. 
Therefore, the magnetic response is diamagnetic. If, in- 
stead of a conventional junction, a 7r-contact is inserted 
in the same sensor, the a(Hd c ) dependence has a local 
minimum at Hd c — and $i > <i> e , providing a param- 
agnetic response. The phase shift a can be obtained from 
Eq. H defining X m = d^./d^e - l. 10 The value of k 2 Q 
allows the estimation of the resolution of the rf measure- 
ment, in this case of the order of a few percent of the flux 
quantum. 

As we discussed above, within the framework of the d- 
wave scenario, extrinsic effects such as faceting can play 
a relevant role in the determination of the properties of 
the junctions and can contribute in particular to causing 
spontaneous currents and/or a paramagnetic effect. In 
the case of a strongly meandering interface, the GB ex- 
hibits a random parallel array of 0- and 7r-contacts. This 
is somehow the analog of the explanation of the PME 
in BiSCCO crystals in the 1-dimensional case of a GB 
line (7r-loops model). An alternative mechanism of the 
PME can be given in terms of the midgap states (MGS) 
and surface properties of rf-wave superconductors (MGS 
model). In this case, the MGS model is valid even for 
flat interfaces of a d-wave superconductor. 

The interplay between the PME and superconductiv- 
ity apart from being an interesting topic itself on d-wave 
induced effects, may be crucial to improve understand- 
ing of transport properties of GBs. It has been demon- 
strated that GB transport properties strongly depend on 
the quality of the substrate^ of the thin film and the 
type of growth. This led to apparently conflicting be- 



(001) YBCO 




FIG. 1: a) Sketch of the investigated biepitaxial grain bound- 
ary where a 45° c-axis tilt accompanies the 45° a-axis tilt; in 
b) a Scanning Electron Microscope image relative to the this 
GB is shown. The elongated grains on the left typical of 
the (103) growth may generate clean interfaces with reduced 
faceting, as evident also by the sketch in c). 



haviors in a wide spectrum of transport regimes. Zig-zag 
Nb-Au-YBCO junctions isolated the effect of facets by 
investigating the presence of spontaneous currents and 
anomalous magnetic patterns^ We intend to reach the 
other limit by investigating a morphology with nomi- 
nal very reduced faceting, and low barrier transparency. 
To this aim we employed biepitaxial junctions, where 
"clean" basal plane GBs can be reproducibly obtained. 

The biepitaxial technique allows the fabrication of vari- 
ous GB JJs by growing different seed layers and using sub- 
strates with different orientations^^ In this experiment 
we have used Ce02 as a seed layer material deposited on 
(110) SrTiOs substrates. Details of the fabrication pro- 
cedure can be found elsewhere^ YBCO grows along the 
[001] direction on Ce02 seed layers, while it grows along 
the [103]/[013] direction on SrTiOs substrates. By us- 
ing Ce02 as a seed layer we were able to induce a 45° 
rotation of the a — b plane of the YBCO with respect 
to the in-plane direction of the SrTiOs substrate, and as 
a consequence the ir contact is formed4& The measure- 
ments shown below are for the case where a 45° c-axis 
tilt accompanies the 45° a-axis tilt (Fig.^i.) This config- 
uration in principle leads to interfaces where effects due 
to faceting can be very reduced, as shown for instance in 
the Scanning Electron Microscope image of Fig. and 
the relative sketch Fig. For biepitaxial junctions in 
the tilt case^ the YBCO growth kinetics and the junc- 
tion interface orientation determine that the long side of 
the [103] grains faces the c-axis counter-electrode, and 
this leads to a more controlled GB (basal plane GB). 
This has been confirmed by cross section Transmission 
Electron Microscope investigations^ 

The sample had originally a typical SQUID geometry 
with a central hole of 50 /im x 50 /im. Rf and SSM 
measurements were performed on a configuration where 
one of junctions was removed and the other was reduced 
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to about 60 fim. 

Results of rf measurements as a function of an ex- 
ternally applied magnetic field at the temperature range 
from T= 4.2 K up to 40 K are presented in Fig. [5] for 
nominally zero field cooled samples (the rest field is be- 
low 0.2 mOe). The sharp minimum at H^—O in Fig. [21 
represents an unusual feature with respect to most GB 
systems, which exhibit a maximum for zero field. The 
minima at a finite magnetic field are originated by the 
redistribution of the magnetic flux into the GB as exten- 
sively discussed in Ref. Il7l 

According to the discussion above, the minimum of a 
at Hdc = is direct evidence of paramagnetic behav- 
ior of this type of GB. Furthermore the paramagnetic 
response was absent after the removal of YBCO from 
the GB region, demonstrating that the effect is only due 
to the GB. In other words we just repeated the same 
measurement on the same sample after removing only a 
narrow region of YBCO along the grain boundary line in 
order to prove that the paramagnetic signal was caused 
by the grain boundary and not by any uncontrolled envi- 
ronmental reason (substrate or sample holder) or possible 
impurities in the YBCO thin film far from the GB line. 
The absence of hysteresis for a^H^-dependence with re- 
spect to external field means that there is no spontaneous 
surface current or flux generated in the GB. This is anal- 
ogous to the situation in rf SQUID with 7r-contact for 
(3 < 1. As a matter of fact, in the case of finite fluctu- 
ations the jumps of $i can occur with a certain proba- 
bility when $ e falls in a hysteresis region— When the 
distribution width of the jump probability is larger than 
the width of the hysteresis, the flux jumps many times 
during the measurements. As a result the apparent time- 
averaged ^(^e) dependence presents a finite slope rather 
than hysteresis even for (3 > 1 (see Fig- - This situation 
has already been observed experimentally— 

A detailed insight into the problem of the param- 
agnetic effect and of spontaneo us cur rents in facetted 
GB's has been given in Refs. I20I2H In this model 
the current density js( x ) ( x is the coordinate in the 
plane of the contact) is a random and alternating func- 
tion of x, jc(x) = (ic)[l + gof(x)] with (f(x)} = 0, 
max[ /(a;)]=l and go — msx\j c {x)/(j c )}. The length 
scale I of g(x) variations is of the order of the grain 
boundary meandering, typically I is in the range of 0.01- 
0.1 /zm— For a conventional current-phase relationship 
j c {x) — (j c ) sin ip(x), the properties of the GBJJ are 
determined by the parameter— 7 = g^l 2 / '8tt 2 A 2 . Here 
Aj = (c$o/167r 2 A(j c )) 1 / 2 is the effective Josephson pen- 
etration depth and A is the London penetration depth. 
Within this model, for 7 < 1 there is no created flux in 
GBJJ and the microjunctions stay in a "excited current- 
less state" — In principle, such junctions can exhibit a 
paramagnetic response. When 7 > 1 there is flux gen- 
erated in the GBJJ. Qualitatively, this requirement is 
similar to the condition (3 > 1 for the rf SQUID, as 
discussed above. Thus the same arguments concerning 
thermal fluctuations can be done for our experiments. If 
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FIG. 2: The phase angle a as a function of the external 
magnetic field applied parallel to c-axes of 65/xm-wide bridge 
across an asymmetric 45° YBa2Cus02, grain boundary. From 
top to bottom the data correspond to T=40, 30, 20, 10, and 
4.2 K. The data are vertically shifted for clarity. 



at <I> e = and T = there is spontaneous flux in the GB, 
one should observe hysteresis of the $i($ e ) dependence 
(see Fig. curve 1). However, at finite temperature 
fluctuations wash out the hysteresis, therefore there is no 
spontaneous current (see Fig. [31 curve 2). The above sce- 
nario relies on a macroscopic approach to the GB, which 
is a random array of parallel and 7r-junctions, i.e. the 
"extrinsic" effect of faceting. 

The other possible scenario is based on microscopic 
arguments, and mainly on mechanisms based on Andreev 
reflection at the surface— iSiSi^Sli The d-wave symmetry 
of the gap leads to the formation of the surface state 
at zero energy (so-called midgap states - MGS2£). This 
MGS is degenerate with respect to the direction along 
the surface of the superconductor (±k y ). Any mechanism 
which is able to split the MGS will lower the energy of the 
system— The splitting of the MGS leads to spontaneous 
currents along the surface and therefore to time reversal 
symmetry breaking. 

Different mechanisms of the MGS splitting were pro- 
posed in the literature (see Ref. |26J and Refs. therein). 
Apart from suggestions based on the presence of an imag- 
inary component of the order parameter (like d ± is, see 
for example Ref. l29h . there are at least two scenarios 
based only on a pure d-wave symmetry of order param- 
eter. 

The first scenario takes into account the Josephson ef- 
fect. For asymmetric 45° GBJJ the energy of the MGS 
can be written as 23 



£mgs(<A k v 



(3) 



where E = A L \A R \D(#)/(2\A L \ + D($)[\A R \ - \A L \}), 
D(d) is the angle-dependent barrier transparency and 
Al, Ah are superconducting energy gaps for the angle 
$ = aicsmky/kF in the left and right superconductor, 
respectively. For asymmetric 45° GBJJ the Josephson 
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FIG. 3: Schematic drawings - time averaged magnetic flux 
inside the GB vs applied magnetic flux (1) - without thermal 
fluctuation, (2) - large thermal fluctuation. 




a) b) 

FIG. 4: Scanning SQUID Microscopy image of a)a 200 x 200 
/tm 2 area, enclosing the tilt GB; the sample was cooled in zero 
magnetic field; b) a 600 x 400 /im 2 area, enclosing the tilt 
GB; the sample was cooled in a small magnetic field. Both 
images refer to the sample measured by rf techniques, and 
the grain boundary junction is indicated by a dotted line. In 
all images, taken at T = 4.2 K, there is no evidence of spon- 
taneous currents. In b) the superconducting nature of both 
electrodes is proven by the presence of Abrikosov vortices, 
which appear elongated in the (103) electrode. 



energy minimum corresponds to the equilibrium phase 
difference across the junction (p = tt/2 which leads to the 
splitting of the MGS. 30 Hence, the spontaneous paramag- 
netic current flows in the surface layer ~ £o- This current 
is screened by the Meissner supercurrent on the scale of A 
and the system pays for the gain in the Josephson energy 
at the expense of the energy of the magnetic field. 

Within the framework of the second scenario the free 
surface of a d-wave superconductor has been found to be 
responsible for the appearance of the paramagnetic effect. 
In this case the paramagnetic quasi-particle current tries 
to compensate the Meissner supercurrent^ As a matter 
of fact the supercurrent causes a redistribution of quasi- 
particles in k-space due to the shift E^ = E^ + Pk ■ v s 
in the quasiparticle excitation energy, where E^ is the 
excitation energy in the absence of a supercurrent, v s is 
the velocity of the supercurrent and pk is momentum of 
the quasiparticle. This is obviously relevant for a (110) 
oriented HTS surface due to the existence of the MGS. 

The key feature for both of the described approaches 
is the presence of a mechanism able to split the MGS and 
produce a gap E$, or in other words to populate the ik y 
MGS unequally. This phenomenon is accompanied by 
a phase transition to a broken time reversal symmetry 
(TRSB) state. The transition temperature Ttrsb can 
be evaluated and in both approaches Ttrsb have been 
estimated to be below 1 K, in agreement with tunneling 
experiments j2i Ttrsb will be apparently the relevant pa- 
rameter also for the PME. However, in analogy with rf 
SQUID properties, above Ttrsb one should observe no 
hysteresis but steep slope of $j($ e ) at <I> e = (again see 
Fig. |3J) . This explains why we observed PME in a wide 
temperature range well above the Ttrsb values but no 
hysteretic behavior. 

Additional information has been provided by SSM in- 



vestigations, which revealed the absence of spontaneous 
currents, as shown in the SSM images of the GBJJ region 
in Fig. E| Fig. Ufa) and Fig. H{b) are SSM images of a 
wide area around the GB "in zero field cooling, and in 
non-zero field cooling respectively, both at T=4.2 K. Fig- 
ure Ufa) confirms in zero field cooling the absence of any 
spontaneous magnetization for this sample. This result 
is different from other measurements on 45° asymmet- 
ric bicrystal 3 ^ and biepitaxial GBJJi2& We attribute this 
difference to the concomitance of reduced faceting along 
the GB and the low barrier transparency of this junc- 
tion, characterized by low values of the critical current 
density jc about 10 2 A/cm 2 . 34 Fig. 0] (b) gives evidence 
of vortices in both the electrodes, and in particular of 
anisotropic vortices in the (103) electrode. As a con- 
sequence, it seems to be not the case that the PME is 
originated by 7r-loops across the GB (7r-loops model, see 
above). Moreover, the observed narrow dip on a vs. Hd c 
characteristics at Hd c — (see Fig. |2} requires 7 « 1. 
Simple estimations show that 7 w 1 corresponds to crit- 
ical current densities across the GBJJ of the order of 
10 5 A/cm 2 , which seems to be unrealistic, given that the 
average critical current density for our junction are of the 
order of 10 2 A/cm 2 . On the other hand the MGS model 
is apparently consistent with our results. 

In conclusion, we have measured the magnetic-field 
response of an asymmetric 45° grain boundary in a 
YBa2Cu30a; thin film. The results of these investiga- 
tions allow us to identify Andreev bound states as the 
cause of the paramagnetic effect, confirming theoretical 
predictions. Andreev bound states have been studied in 
detail theoretically in HTS Josephson junctions and sys- 
tems. However, due to their extreme localization and 
stringent survival conditions, experimental detection has 
been basically confirmed only by one type of measure- 
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ment, zero bias anomalies in tunneling spectral* Our 
method is complementary and direct, and relies on the 
comparison of various types of measurements realized on 
the same samples. The contribution to the existing state 
of knowledge on d-wave order parameter symmetry in 
HTS superconductors is to demonstrate experimentally 
the occurrence of bound states through an innovative ap- 
proach, shedding light on the paramagnetic effect in a 
novel topological configuration. 

A paramagnetic response at low external fields was 



found in a zero-field-cooling regime. We have shown that 
the observed phenomena can be explained by the para- 
magnetic response of MGS surface currents. 
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